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investigation  of  the  electronic  excitation  structure  of  the  transition 


metal  oxides  was  begun  which  included  the  effects  of  correlation  of  the 
excited  electron  with  the  surrounding  medium,  relaxation  of  the  electrons 
around  the  hole  left  by  the  excitation,  and  the  interaction  of  the  excited 
electron  with  the  hole.  Recent  theoretical  and  experimental  work  has  shown 
that  all  these  effects  must  be  included  in  order  to  obtain  a detailed  under- 
standing of  the  excitation  struiyfcjirj*  optLeal^cropertiesa  (ponJt-inuec 
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lattice  dynamics,  and  dielectric  response  of  these  oxides  and  many  other 
materials.  Local  orbitals  instead  of  the  usual  spatially  extended  Bloch 
representation  were  used.  The  local  orbitals  are  especially  useful  in 
describing  the  d-electrons,  which  are  localized  in  position  and  energy. 

The  extensive  codes  necessary  for  the  efficient  calculation  of  multicenter 
integrals  were  constructed  and  assembled  during  the  past  year.  The 
methods  necessary  for  transforming  the  band  states  from  the  Bloch  to  the 
local  orbitals  representation  were  developed,  and  the  methods  for 
diagonalizing  very  large  array  were  investigated. 
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VANADIUM,  AND  THEIR  OXIDES 
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A.  Research  Objectives 


We  began  with  the  following  objectives: 

to  investigate  the  electronic  excitation  structure  of  transition 
metal  oxides 

to  do  the  study  in  the  framework  of  the  local  orbitals  method. 

During  the  past  year  we  have  pursued  the  following  specific  objectives: 
construction  of  calculational  methods 

development  of  codes  for  efficient  calculation  of  individual  multicenter 

two  electron  integrals  using  Slater  Type  Orbitals 

development  of  codes  for  transforming  the  joint  density  of  states 

information  from  a Bloch  representation  to  a local  orbitals 

representation 

development  of  codes  for  reordering  the  integrals  into  groupings 
describing  particle-hole  interactions 

beginning  the  development  of  methods  for  diagonalizing  very  large  arrays 
to  find  the  excitation  spectra  using  an  extension  of  the  methods  of 
degenerate  perturbation  theory. 
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B.  Status  of  Research 


B. 1 Development  of  Calculation  Methods 

As  we  discussed  in  our  original  proposal,  the  excitation  spectrum  is 
contained  in  the  poles  of  the  particle-hole  (PH)  propagator  n_(w).  We  may 
construct  the  PH  propagator  with  the  eigenvectors  C(w)  from  the  matrix  equation, 
represented  here  schematically. 


n(oj)  * £(<d)  ■ A(u>)  £(u>) 


(1) 


with 


n(w)"1  - n°(u)_1  - k(u>) 


(2) 


A(  u>®  ) - 0 

ex 


m 


The  solutions  a>ex  are  the  excitation  energies  which  Induce  poles  in  B (w) . 


The  primary  equation  to  be  solved  is  Eq.  (1)  with  the  structure  presented 


in  Eq.  (2).  The  matrix  11° (w)  ^ has  the  form 


_o  /*-l 
nss’(ta,) 


l Tg*(n1^,n2^)x°(n1k,n2^+q;co)“1  Tg’ (n^.njk+q) 


n^k 


(3) 


with 


T^njk.njk+q)  E e^Cn^+q^e^,  (n^) 

H' 


(A) 


ik-R 


i(k+q) 


J S0  i,cj*  S0  i,  dt,C 


il M. 


s - { j i ;0a  } 


[I 


I 


expansion  of  Che  Bloch  function  and  S is  the  overlap  matrix  between  local 
orbitals.  * is  just  the  inverse  of  the  polarizability  for  the  non- 
interacting particle-hole  pair,  and  it  is  given  by 


(5) 


- f £v*(l“f  £){<*>+  E £ + - E , 
n9k+q  n.k  n,k+q  n.k 


where  En£  and  f £ are  the  one-particle  energy  and  occupation  number  for  Bloch 
state  nk.  It  is  the  inverse  of  Eq.  (3)  that  is  normally  used  in  place  of 
n<w)  in  Eq.  ( 1 ) to  calculate  the  optical  absorption.  The  inclusion  of  K(u>) 
in  Eq.  f 1 > introduces  particle-hole  interactions.  The  first-order  expression 
for  K(u)  in  terms  of  local  orbitals  is 


Kss'  ^aibj ,b' J 'a' i'  Vaia' i' jb' j 'bj 


(6) 


with 


'aibj’.a'i’b'J'  " <(>ai(?>Vi’(r,)  *bj  (r)^b' j ’ (r ' } 


Under  certain  assumptions,  the  one-particle  energy  band  E^  in  Eq.  (5) 
and  the  corresponding  Bloch  coefficients  used  to  construct  Ts(n^k,n0k+q)  in 
Eq.  (4)  may  be  replaced  by  quantities  which  include  correlation  and  relaxation. 
Basically,  these  assumptions  are  valid  when  the  sturcture  of  the  hole  (or 
particle)  is  not  affected  much  by  the  particle-hole  interaction.  We  have 
done  calculations  on  atomic  systems  which  show  that  this  simple  inclusion  of 
relaxation  and  correlation  effects  along  with  including  the  effect  of  the 
particle-hole  interaction  results  in  much  improved  excitation  energies.1 
The  work  of  Kunz  et  £l  on  the  optical  absorption  of  LiF2  suggests  that  in  certain 


cases  such  an  Inclusion  of  separate  effects  is  valid  for  solids.  Certainly 
such  a separation  of  effects  would  be  fruitful  whenever  it  is  suspected  that 
an  excitation  results  in  a well-localized  hole,  as  is  the  case  for  the  d-bands 
in  transition  metals  where  calculations  of  the  dielectric  matrix  have  shown 
that  local  field  effects  are  important.3 

The  solution  of  Eq.  (2)  using  Eqs.  (3-5)  requires  the  diagonalization  of 
a huge  matrix.  However,  if  our  local  orbitals  basis  is  well  chosen,  only  a 
few  terms  from  that  matrix  are  large  for  a given  range  of  excitation  energies. 
To  obtain  a method  which  will  allow  us  to  solve  Eq.  (2)  in  a manageable  sub- 
space while  retaining  the  ability  to  calculate  the  effects  of  terms  not  in  the 
principle  subspace,  we  have  developed  a variation  on  degenerate  perturbation 
theory.  Suppose  the  matrix  jHw)  * is  arranged  as  illustrated  in  Fig.  1.  We 
have  assumed  the  total  matrix  is  of  size  Nm  x Nm  and  have  arranged  the  elements 
along  the  diagonal  into  N mxm  blocks,  II  J ....  TI^  . Then  the  eigenvalue  prob- 
lem becomes 


^ Si  W-f1)*1  \A 


+ i nT1  (A-n_1)~1n“1(A-n  )-1n~jc.  + 

p*v»»l  ~lu  'nJV  ~ 1 1 


(7) 


where  A is  a true  eigenvalue  of  the  total  matrix  and  is  a projection  of  the 
total  eigenvector  on  the  subspace.  The  elements  of  IT  1 could  be  arranged  in 
order  of  Importance  with  respect  to  II  so  that,  for  example,  II  ^ is  larger 
than  II  ^ and  so  on.  In  practice  we  expect  the  series  in  Eq.  ( 7 ) to  be  rapidly 


convergent . 


B. 2 Calculation  Program 


During  the  first  year  of  this  research,  we  have  made  a good  deal  of 
progress  toward  setting  up  our  program  to  calculate  electronic  excitation 
structure.  One  of  the  major  tasks  completed  was  the  construction  of  computer 
codes  to  calculate  general,  multicenter  two-electron  integrals  in  an  STO 
basis.  We  have  also  set  up  the  necessary  programs  to  transform  and  rotate  our 
integrals  between  equivalent  lattice  sites.  These  programs  are  now  running 
and  are  the  only  completely  general  programs  for  calculating  two-electron 
integrals,  actually  without  regard  to  basis  type,  of  which  we  know.  We  need 
these  general  codes  because,  as  is  indicated  in  Eq.  6,  the  two-electron  inte- 
grals do  not  involve  just  matrix  elements  between  one-electron  densities  as 
is  the  case  in  energy  band  calculations.  The  electronic  excitation  problem 
involves  two-electron  density  functions,  and  thus  we  need  individual  two-electron 
matrix  elements.  Fortunately  we  are  helped  by  the  nature  of  the  local  orbitals 
method.  The  method  produces  orbitals  which  are  separated  into  independent  sub- 
spaces over  occupied  and  unoccupied  bands.  Since  v * in  Eq.(5)  involves  only 

— o 

occupied,  unoccupied  pairs  of  Bloch  states,  the  local  orbitals  in  Eqs.  (3)  and 
(5)  occur  only  in  occupied,  unoccupied  pairs,  too.  This  property  greatly 
reduces  the  size  of  our  matrix  problem. 

We  have  also  set  up  the  necessary  computer  codes  to  calculate  the  inverse 
of  the  PH  propagator  given  in  Eq.  (2).  We  use  the  method  outlined  in  the  last 
section.  The  Bloch  sums  involving  the  transformation  matrix  Tg  in  Eq.  (4) 
are  performed  using  the  Lehman-Taut  method.  We  have  generalized  this  method 
to  calculate  efficiently  Bloch  sums  over  the  entire  Brillouin  zone  for  all  the 
combinations  of  bands.  We  can  vary  the  number  of  shells  of  atoms  included  in 
the  calculation,  and  the  central  matrix  which  we  diagonalize  contains  only  the 


most  important  two-center,  two-electron  integrals  for  a given  excitation  energy. 
The  three-  and  four-center  integrals  are  included  as  needed  in  the  perturbation 
expansion. 

We  have  recently  obtained  good,  self-consistent  local  orbitals  energy 
band  results  for  TiO  calculated  by  Jennison  and  Kunz . These  results  were 
obtained  with  the  local  exchange  approximation  rather  than  with  Hartree  Fock 
exchange  and  thus  have  a gap  between  the  p and  d densities  cf  states  in 
agreement  with  the  most  recent  photoemission  measurements.  Evidently  the 
difference  between  Hartree  Fock  and  local  exchange  is  due  to  the  inclusion  of 
relaxation-correlation  effects  in  the  latter's  energy  bands.  Thus  our  optical 
calculations  automatically  include  some  correlation  effects. 

We  are  also  planning  our  calculations  of  the  excitation  structure  in  TiO. 

The  optical  spectra  of  several  of  the  transition  metal  oxides  show  strong 

charge  transfer  structures  just  above  the  band  edge  in  energy.  The  calculated 

d-bands  for  TiO  are  broader  than,  say,  for  NiO.  However,  since  the  final  con- 
++ 

figuration  is  Ti  ,0  , we  expect  charge  transfer  structure  to  exist  in  the 

TiO  spectra.  Ichikawa  £t  al_5  have  also  compared  the  Til.^  X-ray  emission 
spectrum  for  TiO  with  the  M.,  ^ spectrum.  The  initial  states  responsible  for 
the  two  spectra  both  have  p symmetry  and  the  final  states  are  common.  Then 
it  is  expected  from  the  itinerant  electron  model  that  these  spectra  should  show 
similar  profiles.  They  do  not,  and,  furthermore,  the  M.,  ^ spectrum  does  not 
resemble  the  calculated  density  of  states.  These  differences  are  evidently 
due  to  particle-hole  interactions.  We  have  an  approach  to  calculate  optical 
properties  which  includes  particle-hole  interactions  and  is  based  on  local 
orbitals  results  which  give  good  descriptions  of  ground  state  properties. 

Without  such  an  approach,  one  cannot  explain  the  above  localized  excitation 


structures  and,  at  the  same  time,  explain  the  structure  in  the  spectra  due 
to  excitations  of  itinerant  electrons  between  energy  bands. 
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